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Abstract 

In this work a neutrino gas in equilibrium is studied both at T = and at finite temperature. 
Neutrinos are treated as massive Dirac quasi-particles with two generations. We include self- 
interactions among the neutrinos via neutral currents, as well as the interaction with a background 
of matter. To obtain the equilibrium properties we use Wigner function techniques. To account 
for corrections beyond the Hartree approximation, we also introduce correlation functions. We 
prove that, under the quasi-particle approximation, these correlation functions can be expressed as 
products of Wigner functions. We analyze the main properties of the neutrino eigenmodes in the 
medium, such as effective masses and mixing angle. We show that the formulae describing these 
quantities will differ with respect to the case with no self- interactions. 
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I. INTRODUCTION. 



Propagation in dense media is one of the most interesting issues in present neutrino 
physics. The consequences on solar and atmospheric neutrino data, as well as in baseline 
neutrino experiments, are crucial for the understanding of neutrino properties. Moreover, in 
astrophysical and cosmological scenarios, they provide an important feedback which modifies 
the physical evolution of the system under consideration. This, in fact, is the case for 
supernova explosions, where neutrino interactions and oscillations can change the shock 
dynamics 

The nucleosynthesis of heavy elements via r-processes is also largely affected by neutrino 
oscillations, since the proton to neutron ratio can be altered by the oscillations. 

In a similar way, they have been shown to play an important role in establishing the 
cosmic flavor content when the Big Bang nucleosynthesis starts, thus influencing the helium 
production. A crucial ingredient in these two scenarios is the neutrino self- interaction 

mm 

. As it has been remarked, such interactions are non-diagonal in flavor space, and 
give rise to new phenomena in the oscillatory behavior j^, Q| . 

In this work, we analyze in detail the equilibrium properties of a self-interacting neutrino 
gas. Here, we study a system consisting on two generations of neutrinos in equilibrium with 
a matter background and self-interactions taken into account. The formalism, however, can 
be extended in a straightforward way to the case of three generations. 

We assume that chemical equilibrium has been reached, therefore the chemical potential is 
the same for both neutrino species. This is, in fact, the situation for muon and tau neutrinos 
inside a supernova, where both kinds of neutrinos are produced in pairs and, therefore, the 
chemical potentials are zero. It seems also the case in the early Universe, with the present 
values of mass differences and mixing angles, where chemical equilibrium is achieved (at 
east approximately) just before the nucleosynthesis epoch for the three neutrino families 



Neutrinos will be treated in the quasi-particle approximation, i.e. we assume that the field 
corresponding to the mass eigenstates can be described as a superposition of plane waves with 
a modified (with respect to the vacuum) dispersion relation. As we will show, the equilibrium 
features, such as the effective masses and mixing angles, are modified when the interactions 
among neutrinos are considered, due to the non-diagonal nature of the self-term. To this 



2 



end, we use a method based on Wigner Functions, which has been shown to be appropriate 
to describe both the equihbrium and kinematics of many-particle neutrino mixing 
To account for corrections beyond the Hartree approximation, we also introduce correlation 
functions. We prove that, under the quasi-particle approximation, these correlation functions 
can be expressed as products of Wigner functions. As we show, the results obtained using 
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these techniques agree with previous calculations 

This paper is organized as follows. In section II, we derive the equations satisfied by 
the Wigner functions of self-interacting massive neutrinos. In section III we consider the 
Hartree approximation (i.e., when correlations are neglected). In this case, the interactions 
among neutrinos contribute in the same way as the neutral current interactions of neutrinos 
with other particles, such as electrons, protons or neutrons. The addition of correlations 
results in the appearance of non-lineal effects due to the self-interaction. This is described 
in section IV. In section V we analyze the resulting dispersion relations of eigenmodes. Our 
main results are summarized in section VI. The appendix at the end contains the derivation 
of the statistical correlations used in section IV. 



II. EQUATIONS FOR WIGNER FUNCTIONS. 

We consider a neutrino gas consisting on two generations of neutrinos which interact with 
themselves by means of neutral current interactions. We denote the two flavors by (e,/i) . Of 
course, the formalism can also be applied to any two flavors, such as /i and r neutrinos. 

Since we deal with two neutrino species, it is convenient to introduce vectors and matrices 
in flavor (or mass) space. We therefore define the neutrino and antineutrino vector fields: 



, V{x) ^ v\x) lf[x) (1) 



The symbol " on top of a magnitude means that we are dealing with a quantum operator. 
This will be used to distinguish this magnitudes from statistical averages. We also introduce 
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the following matrices in flavor space 

^0 

^0 7^(1 -7^) 

7^(1-7^) 


With the aid of these notations, the Lagrangian terms which account for the neutrino 
self-interactions 0] can be written as 

G p 



ab,/! 



(2) 
(3) 



L{x) = i'{x)i^^dfj){x) — u{x)Mi'{x) 



4^2 



iy{x)il^u{x)iy{x)il^u{x) 



-v{x)K^v{x)v{x)J^fS^{x) - -^^u{x)Q'^u(x)u{x)Afjy{x), 



(4) 



where Gi? is the Fermi constant and M is the neutrino mass matrix. The equation of motion 
then reads as 



iYd^^Hx) - Mu{x) - -^nf'u{x)V{x)n^,u{x) 
V{x)nf'u{x)n^u{x) - -^Af'u{x)v{x)A^u{x) 



4^2 
G p 



AV2 



4^2 



d{x)A^'i}{x)A^j){x) 



2V2 



Vt^V{x)v{x)A^d{x) 



Gf 



2V2 



V{x)n''u{x)A^u{x) = 



(a similar equation can be derived for the adjoint fleld z/(a;) ). 
We now introduce the neutrino Wigner operator 

Ftl{x,p) = {27c)-'l d'y e-^^yd]{x + y/2)i)^{x-y/2) 

From Eq. (jHl) one can easily show that the Hermitian conjugate is given by 

One obtains the following equation for the neutrino Wigner Operator : 

-i[dF{x,p) - 2ipF{x,p)] + 2iMF{x,p) = -(27r)-^^ / dS'd^k e-^''^'^-^'^ 



2V2 

[nF{x,p- k/2)V{y')W{y') + V{y')^v{y')AF{x,p - k/2) 
+ V{y')nu{y')VLF{x, p - k/2) + AF(x, p - k/2)v{y')Au{y')] 



(5) 



(6) 



(7) 



(8) 



where we defined the matrices \1/ = ri + 2A and $ = 2f2 + A. 

We are now interested in introducing statistical averages from the quantum operators 
defined above. These statistical averages are called Wigner functions Q], and are the 
analogous to the distribution functions we need to describe many-particles systems. These 
are, in general, complex functions, and also contain a Lorentz structure. The neutrino 
Wigner functions are defined as: 

i^f F^'{x,p) >= {2ny' J d'ye-'Py < d]{x + y/2)^{x - y/2) > (9) 

Here, the symbol < A > means the average of a given quantum operator A over a basis 
of quantum states which are compatible with the macroscopical knowledge of the system. 
The latter determines a given density matrix operator p . Thus the averaging is performed 
according to 

<A>= Sp{pA} (10) 

where Sp means the trace performed over the quantum basis. Taking into account Eq. 
(|7j) one immediately obtains : 

F\x,p) =j°F{x,p)j° (11) 

which implies that F{x,p) is an Hermitian matrix with respect to generation indices. 

Starting from Eq. (jH)), one can take statistical averages to obtain the equations of motion 
for the Wigner Function, which turns out to be 

[^{d-2ip) + 2iM]F{x,p) = -(27r)-^^ J d^'d^kd^k' 6-"'^''-^'^ 

[ < nF{x,p- k/2)TrmF{y\ k') > + < Tr^F{y\ k')AF{x,p- k/2) > 
+ < TrVLFiy', k')VLF{x,p- k/2) > + < KF{x,p - k/2)TrAF{y', k') >] (12) 

III. NEUTRINO GAS IN THE HARTREE APPROXIMATION 

In this section we investigate the neutrino gas under the assumption of global equilibrium. 
This means that macroscopic quantities must be invariant under space-time translations: 



F{x,p) = F{p) 
5 



(13) 



Moreover, as a first approximation, we neglect the effect of statistical correlations, which 
simply translates into replacing the statistical average of an operator product by the product 
of their averages: 

< Fg'{p)F^t{p') >= Ft^{p)F^f{p') (14) 

Under these conditions, we obtain that the equation of motion for the Wigner Function of 
the neutrinos in equilibrium is [ill 

(7P - M)F{p) = d'kTr[^,{l - ^')F{k)]r{l - ^')F{p), (15) 

In order to obtain a solution to the problem, we simplify our approach to the neutrino 
gas by making the following assumptions: 

1) Neutrinos are assumed to have a small mass (as compared to typical energies). There- 
fore, the neutrino fields consist essentially on left-chirality projections, i.e. we assume that 
the right projections are very small, as compared to the left projections. Moreover, we can 
consider that the neutrino gas contains only neutrinos with negative polarization (negative 
helicity) and antineutrinos with positive polarization (positive helicity). Thus we have, to a 
good approximation, that 

Fip) ~ FLip) ~ F-{p) + F^{p). (16) 

2) Quasi-particle approximation. In the equilibrium gas, neutrinos are in their interacting 
eigenstates, i.e. a particle with definite momentum and polarization has also a definite 
energy. This allows us to treat the neutrinos of the ensemble as free particles with an 
effective mass instead of their own masses in vacuum. In other words, the field can be 
decomposed in plane waves, with a dispersion relation which differs from the one in vacuum. 

3) The mixing angle, which allows us to write a flavor eigenstate as a linear combination 
of effective mass eigenstates, is supposed in principle to depend on all the degrees of freedom 
of the particle, in a similar way to what happens with the mixing angle in matter, when 
neutrinos are propagating trough a matter background, but do not interact among them. 

With all these hypothesis, we can now rewrite the equation of motion as 

(7P - M)Fl{p) = ^/2GFY J d^kTn^FL{k)FL{p), (17) 
or, in a more simplified way, 

(7P - V2GFia - M)Fl{p) = 0, (18) 



where the four-vector a is defined as 

= y (fkTr-f^FL{k). (19) 

In order to obtain an equation which gives the propagation modes of the neutrinos in the 
gas, we multiply the above formula by (7^ — \plGp'^a + M), and neglect the subdominant 
terms proportional to G\. We then arrive to the equation: 

(/ - 2^Gpap - M^)Fl{p) = 0. (20) 

Without any loss of generality, we can take p^ = |"j?|,0,0) , therefore the interaction 
term is ap = a'^p^ — a^\p\. On the other hand, using the known properties of the Wigner 
Functions Fj^ and corresponding to mass eigenstates of neutrinos and antineutrinos, 
respectively, we obtain 

Tr-f'^FLik) = TrYU{k)FLik)U\k) = Tr-ff^FLik) = 

= Tr7'^5^[F--(A:) + F^ (k)], (21) 

a=l,2 

where U is the unitary transformation that transforms from mass eigenstates to flavor states. 
After working out the traces and the integrations, one obtains that the interaction term is 
ap = riyp^ , where riy = rii + n2 = ny^+ riy,^ — — n^^ is the total density of neutrinos 
(with p° ~ \'p'\ ) minus the corresponding antineutrino density. In this way, the equation of 
motion for neutrinos becomes 

{p^ - 2V2GFn,\-p\ - M^)Fl{p) = 0, (22) 

The corresponding equation for antineutrinos (for F^{p) ) can be obtained by replacing the 
— sign into a + sign in the second term of this equation. 

In order to find the effective masses, we have to diagonalize the matrix = ± 
2\/2GF'n^\'p'\. Obviously, this is made by means of the same rotation angle as in vacuum. 
Hence, the effective masses are 

Ml = m\± 2\f2Gpn^\'p\, 

Ml = m\± 2V2GFn^ [pI (23) 

for each generation, where mi and m2 are the vacuum masses, and the signs ± correspond, 
as above, to neutrinos or antineutrinos, respectively. 



Let us now assume that, in addition to the self-interaction among the neutrinos, we have 
an electrically neutral background of matter coniposed by electrons and nucleons. In this 



case the effective masses of neutrinos are 
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li 



(24) 



where 



Ar 



2^GF\'p\ne 
2V2GF\~p\nu 



V2G 



f\ V Fn 



(25) 



being rie the number density of electrons (minus antielectrons) and n„ the number density 
of neutrons (minus antineutrons). 

IV. CORRECTION TO THE HARTREE APPROXIMATION. 

We now want to take into account the effect of the statistical correlations in our treatment. 
For this purpose, we define the two-body correlation function for the neutrino fields as 

Dt^^t{x,x',p,p') ^< F^;>{x,p)F^f{x\p') > -Fff{x,p)F^t{x',p'), (26) 

where subscripts correspond to spin indices. By inserting this definition in the general 
equation of motion Eq. (fT^ . and after some manipulations, it can be written in the form 

[7(9 - 2tp) + 2tM]F{x,p) = -(27r)-^^ / dS'd^kd^k'e-'^^'-y'^ 

2v2 J 

[nTr {D{x,y' ,p - k/2,k')m) + VlF{x,p - k/2)Tr{mF{y\ k')) 

+ KTr{^D{y\x, k',p - k/2)) + Tr {<i>F{y',k')) AF{x,p - k/2) 

+ n {TrnD{y', x, k\p - k/2)) + Tr {nF{y', k')) nF{x,p- k/2) 

+ ATr {D{x,y',p- k/2, k')K) + AF{x,p- k/2)Tr {AF{y', k'))]. (27) 

In Eq. fl27j) . the symbol Tr implies summation over both spin and flavor indexes. 
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As in the previous section, we will assume a situation of global equilibrium, where spatial- 
time invariance is satisfied. For the correlation functions this implies that 

D{x, x , p, p) = D{x — x', p, p') 

When this is applied to the above equation we obtain: 

(7P - M)F{p) = - l')F{p) I d'kTr (7,(1 - 7^)F(A:)) 

+ 1/2 j d^k'd%VtTr {D{k,p- k/2, k')^^ + ATr (^D{k,p- k/2, k')A^ 

+ ATr (^^D{k, k', p + k/2)^ + VlTr (nD{k, k', p + k/2)]^ }, (28) 

where the function D[k,p,p') is the Fourier transformed of the correlation function 

D{k,p,p') = J d'xe-''^D{x,p,p'). (29) 

We will impose the same three hypothesis as in the Hartree approximation, which implies 
having definite mass eigenstates and a mixing angle, which mixes the mass eigenstates to 
produce the fiavor eigenstates of neutrinos. Moreover, under these conditions, we can apply 
the Wick's theorem (see, for example to calculate correlation functions in terms of the 
Wigner Function of the neutrino fields, as derived in the appendix. We obtain: 

DtSik,p,p') = -5\p - p')Ft,j{p' + k/2)F^^{p' - k/2), (30) 

where we are only using the left projections of the neutrino fields in order to construct the 
Wigner Functions, and hence the correlation functions. By inserting the above relation in 
the equation of motion, this one is finally left as 

[7P - V2GFia + V2Gp J d^^FLiqh - M]Fl{p) = 0. (31) 

The four-vector a is defined in the same way as in the previous section. Obviously, the 
third term of this equation provides us with an additional correction to the corresponding 
equation in the Hartree approximation. 



V. PROPAGATION MODES OF THE NEUTRINOS. 



Starting from the latter equation, we can now calculate the propagation modes of the neu- 
trinos in the gas, by performing the following decomposition, consistent with the hypothesis 
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made in the previous sections 

FL{q) = F^{q) + F+{q) = U*{q)F^*{q)U'\q) + U* (q) F+* {q)U*\q) , (32) 

that is, both the corresponding part of neutrinos F^{q) and the corresponding part of an- 
tineutrinos F^{q) of the Wigner Function with flavor indices can be expressed in terms 
of Wigner functions with effective mass indices F£*{q) and F^*{q) for neutrinos and an- 
tineutrinos, respectively (whose mass eigenstates will not be the same as in the Hartree 
approximation), by means of the unitary transformation U*{q) or U*{q), each of them de- 
fined by the corresponding rotation angle. 

On the other hand, we define the following quantities: 

1 r°° 



1 1"°° 

n.,-7r-J d\q\ms'Mq) + c'f2iq)], 

1 r'^ 

n.,, d\q\\q\'cs[h{q) - h{q)l (33) 

where fi{q) and /2(?) are the Fermi statistical distribution functions for each generation, 
corresponding to quasi-particles with well-defined effective masses, and s and c are the sin 
and cos of the rotation angle 9, which relates the eigenstates of effective masses to flavor 
eigenstates. In this way, Ui,^ is the number density of electron neutrinos, n^-^is the number 
density of muon neutrinos and nj,^^ contains interference effects. Analogously, we can deflne 
the number densities for antineutrinos n^^, n^^ and n,yy^. 
Carrying out a straightforward calculation, we arrive to 

J d'qF^iq) = ^iVP^7°, (34) 
where Pl is the left chirality projector and is a matrix deflned (in flavor space) as 

^ Vlf- "Vlf ' '"A, (35, 

\ ni2 n{ii) J ^ \ J ^ \ 2nu S J 

in which n(e) = n,^^ — n,^^ is the electron neutrino (minus electron antineutrino) number 
density, n{ii) — n^,^ — n^^ is the muon neutrino (minus muon antineutrino) number density, 
rii, — n{e)+n{ii) is the total number neutrino density, riu — n^i^— ^v^i ^ — '^(s) ~^(a*) 
is a statistical parameter of asymmetry between the two flavors. It is important to note 
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that the above expressions exactly coincide with those obtained in jl^ by using a totally 
different method. We can finally write the equation of motion in the form 



(7P - V2G'ir7°n^ - v^G'fA^7° - M)Fl{p) = 0. 



(36) 



Obviously, if 5 7^ and ni2 7^ then the interacting mixing angle 6 will not be the same as 
the vacuum mixing angle do. At this point, we can add the contribution from an electrically 
neutral background of protons, neutrons and electrons. This amounts to replacing the matrix 
by 

1 / + S + 2ne — Un 2?2i2 



N* 



2n 



12 



- 5 - rin 



(37) 



In this way, Eq. (jHUj) will now become: 



(7P - V2GFl'^n^ - 2v^Gi.Ar*7° - M)Fl{p) = 0. (38) 

To find the propagation modes, we multiply Eq. (jHHjl by (7J9— v^Gir7°n,^ — 2v^G'i?A^*7°+M) 
, and neglect the terms of order G^., as in the Hartree approximation. Then, one obtains: 

(p2 - V2GFn^P^ - V2GfN*p^ - M^)Fl{p) = 



In order to find the neutrino effective masses we have to diagonalize the matrix 

= M^ + V2GFny + 2V2GfN*p^, 



(39) 



(40) 



and to find the mixing angle in matter we have to obtain the unitary transformation, which 
is determined by the corresponding rotation. After some algebra, we arrive to the following 
expression for the effective masses in the medium: 



Mil = ^ + 2v^G'p (3n, + ne - n„) p' 



=F -A*, 
^ 2 ' 



where 



A* 



(2V2GfP° {ue + 6)- A cos 20o) + (4^2^^^° ni2 — A sin 29o 



1/2 



(41) 



(42) 



is the effective mass difference. In Eq. fl4ip the upper (lower) sign corresponds to M* (M^) 
,where S and A have been defined in Eq. (f^3j) . The mixing angle is given by 

Asm29o-AV2GFP°ni2 



sin 20 



(43) 
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Notice that Eqs. ()41ll43p depend on and 5 which, in turn, have to be evaluated using the 
above equations. In other words, the whole set of equations has to be solved self- consistently. 

Eq. (|39p gives the dispersion relation for neutrinos and antineutrino mass eigenstates, 
which can be written as: 



and provides (as an implicit equation) the energy as a function of the momentum \p\. As 
a first approximation, one can use the fact that, under most situations of interest, neutrinos 
are extremely relativistic particles. Thus, for neutrinos one can replace by \p\. In this 
way, the above dispersion equation can be approximately solved as: 



To obtain the corresponding formulae for antineutrinos we only have to change \p\ to —\p\ 
in the previous equation. 

A consequence of Eq. ()43|1 is that the condition for the MSW resonance is modified with 
respect to the situation where there is not a neutrino background. In fact, the condition for 
the resonance is now: 



This new condition can be of interest if sin2^o — 1, as suggested for both solar and at- 
mospheric neutrino oscillation values. In this case, the MSW resonance might be dominated 
by the neutrino background. In order to investigate this possibility, we need to calculate ni2 
and 6. Using Eqs. and pHjl we arrive to the formulae 



(44) 




(45) 




A [{5 + He) COS 26*0 + 2ni2 sin 29o] 
2V2Gf [{6 + n,f + 4n?2] 



(46) 



ni2 = 



A sin(2^o)^i 
1 + 4:V2GfF2 



S + He 




(47) 



1 + 4:V2GfF2 



where 




(48) 
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The numerator in the integrand of the above equations is a small quantity, due to the 
small mass difference A*(p), therefore it is convenient to expand the numerator using A*(p) 
as a parameter. We then obtain: 

Moreover, in the extremely relativistic limit, we can write 

flip) ^ 



1 + exp [{p - /i) /T] 
In this way, the integrals in Eq. ()48j) can be performed, giving 

Fi = ^ln [l+exp(/i/r)] 

Fi = -^L,2[-exp(/i/r)] (49) 
with Li2 the dilogarithmic function. 

The correction to the MSW condition, as given in Eq. ()4fi|l . is then of the order 

tan/6'o 

ne 

Let us consider the conditions in the Early Universe, just before the nucleosynthesis. We 
then have a temperature T ~ 1 MeV , /i ~ 0.1 MeV and rif, ~ O.lMeV^. By substituting into 
the above expression, one finds that the correction to the MSW condition is only meaningful 
if 

tan 2^0 > 10^^ 

which implies = 7r/4. Such values seem to be disfavored (although not excluded) for 
i^e — ^ i^/iy^r) osciUatious Only if this value is allowed, the neutrino background can 
play a role in establishing the MSW condition. On the other hand, one can check that the 
corrections of the neutrino background to the effective masses, as given by Eq. (jlH), are 
negligible. 

Another possibility consists in the adiabatic Vj. conversion in the presence of a 

neutrino background. The corresponding resonance formula can be found to be 

_ A[^cos2go + 2ni2sin2go] 
^ 2V2G^[52 + 4^y ^ > 

As an example, we consider the proto-neutron star interior with T ~ 40 MeV and yU ~ 100 

MeV. We then find p^ ~ 10^ MeV which is, of course, too high for the neutrinos produced 

in such a scenario. 
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VI. CONCLUSIONS. 



In this paper, we have investigated the equihbrium properties of a system of two genera- 
tions of mixed massive Dirac neutrinos in equihbrium, when self-interactions are taken into 
account. To this end, we have used techniques based on Wigner functions. We assume that 
well-defined quasi-particle states exist for the neutrinos, i.e. the fermion fields for states 
with a definite mass can be expanded in plane-wave states, similarly to the non-interacting 
case, but with a different dispersion relation. 

The equilibrium state is characterized by a single chemical potential /i and a temperature 
T. First, we analyzed the Hartree approximation (when correlations are neglected). In this 
case, self-interactions are diagonal in flavor space and do not modify the mixing angle, 
although they change the effective masses. 

The inclusion of correlations can be done, under the conditions assumed above, using 
the derivations made in appendix. Our results for these corrections agree with previous 
calculations 14, Q], using completely different techniques. These corrections give a non- 
diagonal term in the effective mass matrix. Therefore, in addition to a modification in 
the effective masses of eigenstates, there is a change in the in-medium mixing angle, as 
compared to the Hartree result. Also, the condition for the MSW resonance (when a matter 
background is included) differs from the usual MSW condition. We have shown, however, 
that, for typical astrophysical and cosmological scenarios, the corrections are small, although 
they can be of some interest in cosmology if the i>e-{i>n, Vr) mixing angle is exactly 71/4, as 



still allowed by neutrino oscillation experiments and theoretical models 
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APPENDIX A: CALCULATION OF CORRELATIONS. 

In this appendix we calculate the statistical correlations defined by Eq. (j26p under the 
quasi-particle hypothesis. Let us consider, for simplicity, that we have a neutrino field 
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in equilibrium consisting on one generation of massive neutrinos. Using the corresponding 
unitary transformation to mass eigenstates, the following procedure can be easily generalized 
to more than one generation. The Wigner function for this field is 

Fij{x,p) = (271)-' j Ae-™ < t,(x + 2//2)$,(x-y/2) >, (Al) 

and the correlation functions can be expressed as 

Dijki{x,x',p,p') =< Fij{x,p)Fki{x',p) > -Fij{x,p)Fu{x ,p') (A2) 

We can now evaluate the first term on the right hand using Wick's theorem In this 
way, we arrive to: 

< F,^{x,p)hi{x\p') >= (2^)"'/ dSd^y'e-'^^e-'P'y' (A3) 

< f ,(x + y/2)^i{x - y/2)%i{x' + i/'/2)§,(x' - y'/2) > 
= (27r) y d%d%' e-'P^e-'^'y' [ 

< %j{x + y/2)^i{x - y/2) >< %i{x' + y'/2)^kix' - y'/2) 
- < %j{x + y/2)%i{x' + y'/2) >< - y/2)^k{x' - y'/2) 

+ < f ,(x + y/2)^kix' - y'/2) >< ^,{x - y/2)t,(x' + y'/2) >]. (A4) 

The first term in this expression leads to a product of Wigner functions which is canceled 
by the second term in Eq. ()A2|1 . The second term vanishes. Finally, the third term can be 
expressed as: 

(27r)-« j d^yd^e-'Pye-'P'y' < + y/2)$fc(x' - y' /2) > 

X - < %i{x' + y'/2)^,{x - y/2) > -tSaix - y/2 - x' - y' /2)\ (A5) 

In this formula Sui) is the propagator of the neutrino field. It appears because normal 
ordering of the operators has no been considered. We now impose normal ordering, which 
amounts to neglecting the last term. In this way, we find: 

A,h(x,x',p,p') = -C^^y j d^yd^y'e-'Pye-'P'y' 

< %,{x + y/2)^kix' - y'/2) >< + y'/2)^i{x - y/2) > (A6) 
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By using the equality 

+ z/2)^,{x - z/2) = j d^pe'P'Fji{x,p), (A7) 
and the conditions of equihbrium, which imply that F{x,p) = F{p), we obtain 

D,,kiix-x\p,p') = - j d\5{p-p')e-'''^P'-'^^^^-^'^Fk,{q)Fu{2p' -q). (A8) 
If we take the Fourier transformation of the correlations, we have that 

D,,ki{k,p,p') = {27t)-^J d\x - x')D,,ui{x - x',p,p')e-''("-"'^ (A9) 
and we end up with the following result: 

D,,Mik,p,p') = -5ip-p')Fj,j{p' + k/2)Fuip'-k/2). (AlO) 
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